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Abstract 

We propose a new class of R-parity violating extension of MSSM with type II seesaw mechanism 
for neutrino masses where an unstable gravitino is the dark matter of the Universe. It decays pre- 
dominantly into three leptons final states, thereby providing a natural explanation of the positron 
excess but no antiproton excess in the PAMELA experiment. The model can explain neutrino 
masses without invoking any high scale physics while keeping the pre-existing baryon asymmetry 
of the universe in tact. 
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I. INTRODUCTION 



Recent observation of an excess of positrons in PAMELA experiment [Ij in the energy 
range from 10 to 100 GeV has confirmed similar observations several years ago by the 
HEAT {2! and AMS experiments. A great deal of discussion is currently under way to 
understand this excess and its possible implications for physics beyond the standard model. 
While it is quite possible that this excess is of pure astrophysical origin [3], there is hope that 
this may be coming from the dark matter in our galaxy. Many models have been proposed 
to explain this excess in terms of different kinds of dark matter. 

Interpretation of the observations in terms of dark matter raises several issues: 

• how to explain the lack of any excess in the hadrons? 

• how does one get an adequate enough positron production rate to explain the excess? 

These issues have been discussed in two broad classes of models for dark matter: (i) stable 
dark matter pair annihilation or (ii) decaying dark matter [6|. 

In the first category of models with thermal production of dark matter relic density 
in the early universe, the dark matter annihilation is constrained by the observed relic 
density today. If one uses the same value for the annihilation cross section, an additional 
enhancement is required to understand the observed positron excess for popular dark matter 
density profiles, e.g. Navarro-Frenk- White (NEW) [7| type. There have been two ways to 
address this problem: (i) by increasing the local density of dark matter and/or (ii) by 
introducing new hitherto unobserved light particles [sl that can give rise to the so-called 
Sommerfeld enhancement of cross sections at lower particle energies. 

On the other hand, in the alternative scenario involving decaying dark matter, there is no 
correlation between the cross section that generates the relic density and the magnitude of 
positron excess since the latter involves the decay rate which is independent of the physics 
involved in relic density generation. Therefore the second problem does not arise in this 
class of models. Eor this reason, we focus on a decaying dark matter model in this paper. 

The issue of lack of hadrons can potentially be a problem in both classes of dark matter 
models and requires further model building with specific dark matter properties. 

The decaying dark matter we consider is an unstable, long lived gravitino which can arise 
naturally in supergravity models if gravitino is the lightest super symmetric particle and if 
the model violates R-parity. The long lifetime of gravitino required for it to be a viable 
dark matter is automatically satisfied since the gravitino decay involves a combination of 
the very weak gravitational interactions as well as weak R-parity breaking Absence of 
lepton number violation in any observed process to date justifies the second assumption. 



Let us start by discussing the familiar R-parity violating models [10|, i.e. the minimal 
super symmetric standard model (MSSM) with R-parity violating terms in the superpotential 
of type LLe'^, QLdP, u'^d'^dP and LHu- In these models, one generally expects both leptons 
as well as hadrons in the final states of gravitino decay. Models of this type have been 
considered in ref. [ll[. If however one kept only the LLe'^ term and drop all the others, 
then, the predominant decay mode of the gravitino will only be to leptons as required 
to understand the PAMELA data. However, one problem with this scenario is that if the 
strength of this coupling (usually denoted by A), is larger than 10~^, this will erase the baryon 
asymmetry of the universe. While the strength of this magnitude is also what is required for 
understanding the PAMELA observations, it is too small to explain neutrino masses via loop 
corrections, without assuming other physics beyond MSSM, e.g. grand unification. Thus if 
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we want a minimal bottom-up approach to understand the PAMELA observation, while at 
the same time keeping the baryon asymmetry of the universe untouched and an explanation 
of neutrino masses with only TeV scale physics, one must seek alternative models. This is 
what we do in this paper. 

We propose a new class of R-parity violating interactions that can arise in extensions of 
MSSM which does three things using only TeV scale physics: (i) it explains small neutrino 
masses and mixings via the type II seesaw mechanism [12| ; (ii) it keeps the baryon asymmetry 
of the universe untouched and (iii) it is able to explain the leptophilic nature of the PAMELA 
observations, as a result of gravitino dark matter decay. We also point out that for a different 
choice of parameters of the model, consistent with our other requirements, it can explain 
the recent FERMI observations. 

We also point out a novel feature of any decaying dark matter model which decay to 
photons (even as a subdominant decay mode) that they can be used to map out the dark 
matter density in the galaxy. 

This paper is organized as follows: in sec. II, we present the R-parity violation model 
and discuss some of its implications; in sec. Ill, we consider the gravitino as the dark matter 
and address its decays to show that the model naturally predicts only three leptons as the 
dominant decay final states; in sec. IV, we present our fit to the PAMELA data; and we 
conclude in sec. V. 

II. NEW R-PARITY VIOLATING MODEL 

We extend MSSM by adding a pair of SU{2)i triplets A, A with hypercharge Y = ±2. 
The A field couples to leptons generating neutrino masses when the triplets acquire small 
vacuum expectation value (vev). We include only the new R-parity violating interaction that 
involves the A field and no others. In the appendix, we show that this model is radiatively 
stable. 

The superpotential in our model consists of three parts: W = Wmssm + + 5W^, 
where Wmssm is the familiar R-parity conserving MSSM superpotential; 5W includes the 
new R-parity conserving terms that involve the A and A fields and the last term is the 
term. More explicitly, 

Wmssm = XuQ'^iT2H^u'' + XdQ^iT2Hdd'' + A^L^zrsi/de' + fiHuHa , (1) 

with usual soft terms. 

6W = fL^iT^AL + edHjtT2AHd + euH^ir^AHu + /lATr (AA) . (2) 

The new soft supersymmetry (SUSY) breaking terms associated with 6W are 

6Cs = fAL^ir^AL + edAHjiT^AH^ + e^AHliT^AH^ + h^Ti (AA) + h.c. . (3) 

Since the terms in 5W will be responsible for small neutrino masses, we expect the dimen- 
sionless coupling parameters in 5W to be very small. Roughly speaking, the corresponding 
soft SUSY breaking parameters are of the order euA ~ ^u^s, ^dA ~ ^dMs, with Ms being 
the SUSY breaking scale. 

Let us now discuss the R-parity violating interactions in our model. Note that Hu,d have 
no lepton number whereas A has L = 2; the above superpotential therefore breaks lepton 
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number by two units but does not break R-parity We now add the following interaction 
which violates R-parity: 

SW^^ = aAHdL. (4) 

The associated soft breaking terms is given by: 

Cf^, = pLAHd + h.c. (5) 

Thus the two R-parity breaking parameters in our model are given by a and p. We forbid the 
R-parity breaking terms of MSSM, i.e. XLLe'^ + X'QLd!^ + \"u'^d'^d'^ + p'LHu from appearing 
in the superpotential. The question of radiative stability of this choice of R-parity breaking 
is discussed in the Appendix. We will call this type II R-Parity breaking (RPBII) as opposed 
to the MSSM with R-Parity breaking which we will call type I R-Parity breaking (RPBI). 
It is easy to see that unlike the RPBI models, RPBII models are quite well hidden in low 
energy particle physics processes, even for R-Parity breaking couplings of 0{1). 

We have found a symmetry which will forbid the LLe^ and the u'^d'^dP terms of MSSM 
while allowing the ALHd term. As far as the QLd'^ term is concerned, once it is set to zero, 
a very small value for this coupling is induced in higher orders after A field takes a vev and 
it does not affect our results. 

As far as gauge coupling unification is concerned, if there is no new physics in the theory 
above the TeV scale, the couplings do not unify. However, with extra intermediate scale 
particles, one can restore unification of couplings. We do not address this issue in this paper. 



A. Baryon asymmetry and neutrino mass with R-parity violation 

As noted in the introduction, the leptophilic nature of the PAMELA data could also be 
explained by keeping LLe"^ terms of MSSM but the coupling strength of this interaction A 
must be below 10~^ for both baryogenesis protection 13|] as well as PAMELA explanation. 
They will then lead to one loop neutrino masses of order }?mi/ (IGtt^) ~ 10~^ eV which are 
much too small to explain observations. 

If one uses grand unified models such as SO(IO) to generate neutrino masses via higher 
dimensional operators such as 16ml6ml6Hl6H, then in this model the R-parity breaking 
interactions can come from A16ml6ml6ml6// type terms after B-L symmetry is broken by 
the vev of v'^ field in 1%h- One could then make this couplings to be of the right order, i.e. 
\{y'^)/Mpi ~ 10~^. The problem here is that it also generates u'^d'^d'^ type R-P violating 
terms with similar strength and the model runs into conflict with proton decay. Similar 
situation happen in the SU(5) model i.e. the interaction lOmSmSm that generates LLe'^ type 
R-P violating term also generates u'^d'^d'^ term and therefore will have the same proton decay 
problem. 

The only model where only LLe'^ and QLd'^ terms can be generated without generating 
the u'^d'^d^ term is the model based on the gauge group SU{2)l x SU{2)ji x SU (4)c where 
they arise from separate higher dimensional terms- the first two from FFF'^F'^ and the last 
one from p'^p'^p'^p'^ where F = (2,1,4) and F'^ = (1,2,4) are the representations of the 
gauge group that contain fermions. 

Our model, on the other hand, is a type II seesaw model with all particle masses in 
the TeV range. The constraints on the parameters of such models such that they do not 



erase any pre-existing baryon asymmetry have been discussed in |lj]. The discussion for 
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our model is similar to this paper. A combination of the /, hu,d,i with either eu,d or a will 
lead to violation of lepton number, L. Therefore if the strengths of the L-violating couplings 
are such that the processes caused by them are in equilibrium in the early universe, it will 
erase any pre-exisitng baryon asymmetry of the universe. From [3] we see that for neutrino 
masses in the sub-eV range, the constraints on the triplet vev depend on the mass of the 
A particle and for few hundred GeV mass of A, vt < 10""^ GeV. We can therefore choose 
the parameters of the model such that any pre-existing baryon asymmetry is not erased by 
them. 



B. Symmetry breaking and triplet vevs 

After electroweak symmetry is broken by the doublet Higgs fields {Hu^) = v^^ 7^ 0, A 
acquires an induced vev, f a- The magnitude of triplet vev's can be estimated by minimizing 
the potential of the theory to be: 

(A) = ( ) ' = ~~2 {^dAvl + e*/iAf^ + ed^*VuVd) , 

(^) = Q ^ , Vt = i^uAvl + e>Afd + ^ufJ'*VuVd) , (6) 

where mo is the typical SUSY breaking mass and all other parameters are defined in the 
text. Rough order of magnitude of vt is ~ eu^v"^^/ Ms- If vt < MeV, < 10~^. The 
triplet vev of an MeV is an input into our model. It could arise from tadpole terms as in 
the usual type II seesaw models i.e. by minimization of the terms M^A'^A + \m^i2^HuHu 
by choosing A of order 10"'^. 



C. Neutrino mass and constraints from lepton flavor violation 

Neutrino masse in our model are given by type II seesaw 1^ . 



= 2/t;T ~ O.leV, (7) 

then implies that if vt < MeV, / > 10^''. For vt 7^ 0, there is a — if^ of magnitude avT 
induced which via a seesaw-like formula give an additional contribution to neutrino mass 
~ the choice of parameters in our model, this contribution to neutrino mass is 

MsuSY ' 

< 10~^ eV and is thus negligible. 

We now turn to the lepton flavor violation constraints e.g. — >■ 3e on the couplings /. In 
a generic triplet models of this type, /x — > 3e decay can arise via the exchange of the doubly 
charged component of the triplet A. Present upper limits on this process restrict the values 
of the triplet coupling as follows: /11/12 < 10"^. We will keep the / coupling in the range 
10~^ < / < 10"'^. We will see later that for natural values of parameters of the model, / 
will be closer to 10~^ — 10~^. 

It is also worth pointing out that in this model after standard model symmetry breaking, 
the sneutrino field will have a vev; but this vev has a magnitude {y) ~ avT ~ 10~^^ GeV 
which is much too small to affect any low energy leptonic physics e.g. neutrino mass. 

Furthermore, in our model, we do not have a massless Majoron because we have explicit 
lepton number violation by the terms eu,d- Also since the triplet masses are in the 10^ GeV 
range, there is no new contribution to the Z-width. 
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III. GRAVITINO DECAY AND LIFETIME 



To discuss dark matter and its application to PAMELA, we choose gravitino as the 
hghtest supersymmetric particle (LSP) and hence the candidate for dark matter. If the 
theory conserved R-parity, the gravitino would have been stable. However, since our model 
has R-parity breaking, it is unstable and as we will see below, it will have a long lifetime so 
that it can be a viable dark matter of the universe. To estimate the gravitino life time, let 
us look at its various decay modes and the total width. 

The gravitino can have both two and three body decays. First we will show that the 
two body decay rate is very small for the choice of parameters in the model, making the 
three body decays dominant. To see this, we note that in this model at the tree level after 
electroweak symmetry breaking, A — £ appears, while the other mixings such as W~ — i~ 
and Z — V are severely suppressed by C(e/a) -C l/(167r^) comparing with A — £. The A — £ 
mixing gets contributions from both supersymmetric as well as SUSY breaking parts and is 
given by 

UeA - — 2 r~ ■ (8) 

m~ — 



G wwvx. 




FIG. 1: The Feynman diagram for gravitino decays to three leptons through a triplet Higgs. 



The three body tree level decay of the gravitino arises from the Feynman diagram at 
Fig. [Hand the decay width for G — ^ e'e^Ui is given by 



1927r3 SMp^i 



where x = m\/w?~ and the function F{x) is given as F[x) = {2x — l)(30a;^ — 66x + 37)/12 + 
{x — 1)^(52; — 1) ln[{x — l)/x]. For > mg, the maximum value of F{x) is about 0.04. 

For gravitino as the dark matter, we choose its lifetime to be ~ 10^^ sec, which corresponds 
to the three lepton decay width F ~ 10^^*^ — 10^^^ GeV. For G mass to be ~ 300 GeV, and 
mA > m^, this requires |f/eA/| ~ 10"^. 

Let us now turn to the two body decays of type G —>■ u + 'j,£ + W,!/ + Z etc. These 
will arise from R-parity violating mixings between 7 — z/, W — i, Z — v etc. In our model, 
these mixings are either absent or severely suppressed by Oieja) ^ l/(167r^) at tree level. 
At loop level, they appear in similar way as shown in Fig. [2l To give a typical estimate, we 
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FIG. 2: The one-loop contributions to the R-parity violating mixings between 7 — 1^, W — I, Z — v. 



find for 7 — z/ 



The two body decay width of the gravitino is given by: 

i two body i287rM|,, • 
The ratio of two to three body decay rates can then be given by: 



(11) 



two body 



J- three body 



(12) 



For ^ mA — '^17', there is a parameter range where the three body decay dominates. 
As an example, we choose / ~ 10~^ and a ~ 10~^, we find that If/jy-^l ~ 10~^^ — 10~^^ or so. 
To get |/?7eA| — 10^^ then implies that we have to choose — to about 10% accuracy. 



IV. DIFFUSION AND PAMELA OBSERVATIONS 

In this section we study the positron signals from the gravitino decay. In this model, 
the gravitino dominantly decays to three body leptonic states that involve electrons, muons 
as well as tau's with triplet Higgs mediated. To see this note that the parameter a and 
p in Eq. H] and [5] have fiavor index. For simplicity we could choose this to be along the 
electron fiavor direction. For the case of degenerate neutrinos, then, the gravitino will have 
the following final states: 

G e+(z/ee", i/^/i", i/^r") (13) 

In our fit we take all these modes into account. The positron fiux from the decay of dark 
matter in the halo can be obtained by solving the steady propagation equation 15, 16| : 



V ■ (K(E,f)V/e+) + ^ (6(i?,f)/e+) +g(i?,f) = , (14) 

where /e+ is the number density of positron per unit energy, K{E, x) is the diffusion 
coefficient, b{E,x) is the rate of energy loss and is assumed to be spatially constant as 
b{E) 10~^^(£'/lGeV)^sec~^. The source term Q{E,x) is given by 

M^L^^, (15) 

mgTg dE 
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combining the distribution profile of dark matter and the energy spectrum of positron from 
gravitino decay. The solution of the transport equation at the Solar system can be expressed 



by the convention [11 



Emax dN 4 

dE'G{E,E')^ 



where the Green's function is well approximately given by 

10^6 



G{E,E') 



E^ 



e{E' - E) cm"^sec , 



where the energy is in units of GeV and we adopt parameters a = —1.0203, b 



0.70 with assuming the NFW profile and the MED diffusion model [H 



The positron flux from gravitino decay can then be obtained from 



e+ 



Air 



feAE) 



dE'G{E, E 



dE' 



(16) 
(17) 

-1.4493,5 = 
(18) 




FIG. 3: Fraction of positron flux as a function of energy with the contributions from gravitino 
decay. 



To get the positron flux fraction, one needs to know the background fluxes of primary 
and secondary electrons and secondary positrons. We use the parametrizations obtained in 
0,1 Il7l | with the fluxes in units of (GeV ^ cm ^ sec ^ sr ^): 

16£'~^-^ 



<^lt''{E) 



1 + 11^0-9 + 3.2^2-15 

0.7^0-7 



1 + 110^1-5 + 600^2-9 + 580^4-2 ' 

^TAE) = . , , , (19) 



1 + 650^2-3 + 1500^4.2 
where E is expressed in units of GeV. The fraction of positron flux is then given by 



_e2 

e+ e+ e~ e 



(20) 
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where k = 0.88 is a free parameter used to match the data when no primary source of e"' 

flux 



For simphcity, we consider the degenerate neutrino mass spectrum, in which case the 
triplet Higgs mainly couples to lepton pairs with same flavor. We analytically calculate the 
electron, positron energy spectrums for G ue^e~ as well as their boosted spectrums from 
fi'^ fi~{T~^T~) cascade decays. To fit the PAMELA'S data, as an example, we take mg = 350 
GeV, rriA = 700 GeV, \fUj^\ = 2.5 x 10~^, therefore the lifetime of gravitino is about 
2.1 X 10^^ sec and the fitting result is shown in Fig. [31 



Let us now briefly comment on the recent Fermi-LAT observations [19|. The Fermi 
observation shows a less pronounced e~ + e+ excess above the background in the ATIC 
energy range of above 100 GeV. We believe that one can get a fit to the Fermi data if we use 
normal neutrino mass hierarchy and suppress the RPV coupling involving electron superfield. 
The point is that in this case, the final electrons arise mostly from the decay of the final state 
muon and tau coupling to the A field and therefore have a much softer spectrum. Without 
suppressing the RPV coupling involving electron superfield, one can choose the parameters 
as = 3 TeV, = 2.9 TeV, \Uj^\ = 3.5 x 10"^ therefore = 0.42 x lO^^ sec. In this 
case the gravitino decays to leptons and on-shell A's which will mainly produce muons and 
tauons (then cascade decay to electrons) by choosing the normal neutrino mass hierarchy. 
The positron fraction and electron-positron spectrum with the contributions from gravitino 
decay is given in Fig. HI where k = 0.72 is used to normalize the background. We see that 
for this choice of parameters, both PAMELA and current FERMI data can be fitted well. 




PAMELA 
background 
gravitino decay 




Energy [GeV] 



Energy [GeV] 



FIG. 4: (left) Fraction of positron flux as a function of energy and (right) the electron-positron 
spectrum fit after Fermi-LAT data with the contributions from gravitino decay. 



V. CONCLUSION 

In conclusion, in this brief note we have proposed a particle physics interpretation of 
PAMELA positron excess in terms of a new class of R-parity violating models (RPBII) 
which is related to the neutrino mass via type II seesaw mechanism. This provides a natural 
explanation of why the PAMELA excess is only in positrons and and not in hadrons. For 
the case of normal neutrino mass hierarchy and a different choice of gravitino and Delta 
field masses, the model can describe also the new FERMI data. This model also provides 
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an explanation of both small neutrino masses without invoking any physics above the TeV 
scale while at the same time making sure that any pre-existing baryon asymmetry of the 
universe is not erased by the R-parity violating interactions. This class of R-parity breaking 
models turns out to remain very well hidden from low energy experimental probes unlike 
the MSSM R-Parity breaking models. 

An interesting possibility for decaying dark matter with a photon in the final state is 
the opportunity to measure the dark matter density in the halo using the intensity of the 
gamma rays from different directions in satellites. We will elaborate on this idea in a future 
pubhcation. 
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Appendix: Radiative stability of the vacuum 

In this appendix, we address the issue of whether the smallness of e or eA which char- 
acterize the lepton number violating terms are stable under radiative corrections. /^From 
the discussion of neutrino mass, we know from Eq. ([7]) that e ~ 10~^^/ /, while in order to 
explain PAMELA anomaly, we demand a ^ 10~^//. To suppress gravitino two-body decay, 
we need to stabilize the hierarchy to ensure e ^ a/(167r^). 

For this purpose, note first that radiative correction to e is safe because it is a super- 
potential parameter and is therefore only multiplicatively renormalized due to the non- 
renormalization theorem of supersymmetry. However the soft term is not multiplicatively 
renormalized and there is a contribution which is not proportional to itself, as shown in the 
following figure. 

<H ,> 

d 

I 

I \ 



L -^^^^ E^- 
I 

X 

<H ,> 

d 

FIG. 5: The radiative contribution to ff^Aff^. 

We show what the renormalization of would be by a symmetry argument. First, we 
are considering H^I\H^ soft term, so the radiative correction must be proportional to one 
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of the soft parameters. Second, HaAHd violates lepton number by 2 units, so the radiative 
correction, if not proportional to itself, must include two AL = 1 processes, i.e. oc a^. 

Third, we can extend and restore the Peccei-Quinn (PQ) symmetry by assigning charges 
to fields and spurion's parameters. In MSSM, the PQ charges are = Ye<^ = Yq = Yjjc = 
Yijc = 1 and Y^^ = Yh^ = —2. The PQ symmetry is restored if the parameters also carry 
charges as Y^ = Yi, = 4. Now we extend it and set the PQ charges for triplet Higgs as 
Y^ = Y^ = X. From the fLAL term, we need assign Yf = —{x + 2). In a similar way, 
Y^^ = —2x, Ya = Yp = 1 — x,Y^ = = 4 — X demanded by the terms /iATr(AA), LAHd 
and HdAHd- From the second comment, to renormalize eA whose charge is 4 — x, we must 
need to have two a's, which totally contribute 2(1 — x) charge. Therefore it is still short of 
x + 2, which is exactly the charge of /* without other choices. 

Therefore, one can conclude that the additive renormalization of the parameter must 
take the following form 

JlA^^A^aT, (Al) 
\Vu,kJ 167r^ 

which is 10"^*^// <^ ^a/vwu ~ e- So the smallness of is stable under radiative corrections. 
We also note that the new R-parity breaking superpotential 

Wfi = aLAHd (A2) 

introduces soft R-parity breaking terms of usual MSSM variety from Fig. [61 This term will 
generate the usual R-parity violating MSSM terms through radiative corrections (except the 
A" term). Their strengths are, however, very weak and do not lead to any observable effects. 




FIG. 6: This diagram generates the L — mixing as well as other soft R-parity breaking terms 
of MSSM. 
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